ON FACTORIZATIONS OF SMOOTH NONNEGATIVE 
MATRIX- VALUES FUNCTIONS AND ON SMOOTH 
FUNCTIONS WITH VALUES IN POLYHEDRA 



N.V. KRYLOV 

Abstract. We discuss the possibility to represent smooth nonnegative 
matrix-valued functions as finite linear combinations of fixed matrices 
with positive real-valued coefficients whose square roots are Lipschitz 
continuous. This issue is reduced to a similar problem for smooth func- 
tions with values in a polyhedron. 



1. Motivation 

One of the main goals of the article is to understand what kind of optimal 
control problems of diffusion processes is covered by the results of [3] and [7J , 
where the processes are given by ltd equations in a "special" form, such that 
in the corresponding Bellman equation the second order part is represented 
as the sum of second-order derivatives with respect to fixed vectors (indepen- 
dent of the control parameter) times squares of real- valued functions that are 
Lipschitz continuous with respect to the space variables. Roughly speaking 
the answer is that all control problems with twice continuously differentiable 
diffusion matrices fall into the scheme of [3] and [7J whenever property (A) 
holds: these matrices for all values of control and time and space variables 
belong to a fixed polyhedron in the set of symmetric nonnegative matrices. 
In the author's opinion the control problems with property (A) are the only 
ones which admit finite-difference approximations with monotone schemes 
based on scaling of a fixed mesh. 

For functions w(z) given in a Euclidean space and vectors £ in that space 
set 

i i,j 
In many situations one needs to represent a d x d nonnegative symmetric 
matrix u as the square of a matrix or more generally as the product vv*, 
where v is not necessarily a square matrix. If u = = vv* and v = (y ik ) 
and for each k we introduce the vector v k = (v tk ) £ M. d , then for any smooth 
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fix) given on M. d and the operator 

/./ : Y. ( L1 ) 

hi 

we have 

L f =/ J f(v k )(v h )- I 1 - 2 ) 

k 

In fact, as is easy to see having (jl.2p for all smooth / is equivalent to the 
validity of the formula u = vv* . 

There are very many v such that u = vv* and then a few questions arise: 

(i) if u is a measurable function of a parameter, can one find a measurable 

v? 

(ii) if u is smooth, can one find a Lipschitz continuous v? 

The answer to the first question is easy and positive. Indeed, one can 
take 

poo 

v = c r^ 2 (e~ ut -l)dt, 
Jo 

where c is an appropriate constant. This defines v as the square root of 
u. Since long ago it is known that the square root of a twice differentiable 
nonnegative matrix- valued function is Lipschitz continuous (see [3], |10|). 
This result was used in the investigation of solvability of degenerate ellip- 
tic and parabolic second-order equations by using probabilistic or classical 
approaches. 

However, there are applications in which formula (II. 2D is not very conve- 
nient. One of these applications is related to finite-difference approximations 
of solutions to elliptic and parabolic equations with variable coefficients u 13 . 
Formula (II. 2D suggests replacing f( v k^ v k^ with the second-order difference 
along vectors v k and if v k vary, it may be impossible to find a reasonable 
mesh on which the approximation operator would make sense. This prob- 
lem does not arise if v k = Ckj k , where ^ k are constant vectors and are 
real-valued functions, because then 

Lf = } ] crlf(^h )h h) (1.3) 

k 

and one can concentrate on meshes that are obtained by contracting 

n klk ■ n k = 0, ±1,...}. 

k 

According to Remark 2.1 of [3] considering operators L in form (jl.3)l is 
rather realistic from the point of view of numerical approximations. It turns 
out that if we fix a finite subset B C M d , such that Spanl? = M. d , and if L 
from (jl.ip admits a finite-difference approximation 

L h f(0) = Y J Ph(y)f(.x + hy)^Lf(0) as h [ 0, V/eC 2 

yeB 
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and Lh are monotone, that is Lhf(0) > whenever f(x) > /(0) on M. d , then 
automatically L is written in the form (|1.3p with some 7 fe E B. 

Now the question is: If u = u(x), under which assumptions can one find 
constant vectors 7fc's and functions cf\{x) in order for (II, 3p to hold? Perhaps, 
Motzkin and Wasow (see [9]) were the first to address this question in the 
framework of finite-difference approximation. They proved (see also Lemma 
17.13 of [5]) that if we denote by S[X, A] the (closed) set of positive d x d 
matrices with eigenvalues lying in the interval [A, A], where < A < A, 
then there exist a finite set of unit vectors 7 1 ,...,7 m E M d2 and numbers 
< A* < A*, such that any u E S[\, A] can be written in the form 

m 

u ij = J2^1 ik J jk , (1-4) 

k=l 

where the numbers satisfy the inequalities A* < (3^ < A*. In that case 
(Ol) holds with of = f3 k . This fact was used in the development of the 
theory of fully nonlinear second-order elliptic and parabolic equations. 

One can give a quite easy explanation of this result. If we take any Ai < A 
and Ai > A, the set S[\\, Ai] will contain an open polyhedron P° containing 
S[X, A]. Each point of a polyhedron is represented as a convex combination 
of its vertices and one easily obtains (jl.4p . for instance, as in the proof of 
Lemma 5.5.4 of |6]. 

With a little more effort one can get more convenient representations. We 
claim that given any open bounded polyhedron P° in a Euclidean space E 
of points y with vertices, say yi,....,y n , there exist infinitely differentiable 
functions Pk{y) > 0, k = 1, ...,n, such that for any y E P° 

y = ^2vk(y)yk, ^2pk(y) = l- (1.5) 

k k 

This fact is proved by induction on the dimension of P°. First, without losing 
generality one may assume that the volume of P° in E is strictly greater than 
zero. Then, assume that the fact is true for any face of P° and then solve 
Laplace's equation Ap^ = in P° with boundary condition pk(y) = Pk(y) 
on dP°, where pk is the weight of the vertex in the representation of 
y E dP°, which is supposed to hold by the induction hypothesis. Of course, 
if y E dP° and y^ do not belong to the same face, we set pk(y) = 0. Then 
by the well-known properties of harmonic functions p^ > in P°, they are 
infinitely differentiable in P° and since 

^^Zvk{y)yk = 0, in P° and ^Pk{y)yk = y on dP° 

k k 

and Ay = 0, by uniqueness we have the first relation in (jl.5l) . The second 
one is obtained similarly from the fact that it holds on dP° and Al = 0. 



-1 



N.V. KRYLOV 



After having proved the claim we return to the original P° and write for 
any u € P° 

u = s ^p k (u)u k , y^ y p k (u) = 1, (1.6) 

k k 

where Uk € P° C S[X\, Ai], p\~ are infinitely differentiable in P°, in particu- 
lar, in S[X, A], pfc > in P°, in particular, they are bounded away from zero 
on the closed set 5 [A, A]. Now to obtain (|1.4jl from (|1.6j) it only remains to 
recall that if z = 1, d, are unit eigenvectors of u k with eigenvalues 
then Ai < fiki < Ai and 

i 

The above construction of Pk(u) has a very substantial advantage over 
the original one in [9] (or [5] and [6]). Namely, it is seen that if u = u(x) is 
a smooth function of a parameter x, then in the representation 

u ( x ) = ^2pk{u{x))u k (1.7) 

or in the implied representation (jl.3p the functions pfc(u(x)), p fc (u(x)), and 
(jfe(x) are as smooth as u(x) is. 

We see that from the point of view of the possibility of applying numerical 
approximations to uniformly nondegenerate equations the situation looks 
quite promising. For degenerate equations and fully nonlinear equations the 
situation is much more complex. In this case we again may try to prove 
(|1.6j) with pk such that p k (u) is Lipschitz continuous in u. However, this is 
impossible even if d = 1 and S[\ A] = [0, 1]. In this case, naturally u\ = 0, 

1/2 

U2 = 1, and P2(u) = u, so that p 2 (u) is not Lipschitz continuous. 

On the other hand, in numerical approximation or probabilistic approach 

1/2 

one needs p k (u(x)) to be Lipschitz continuous function of x, rather than 

1/2 

p k (u) to be Lipschitz continuous function of u. This slight difference makes 
the problem solvable in some cases. For instance, in the above case that 
d = 1 it is known that for any nonnegative twice continuously differentiable 
function u(x) its square root u l l 2 (x) is Lipschitz continuous. 

Another example is given by the functions with values in the set of the 
so-called diagonally dominant nonnegative symmetric matrices, which are 
quite popular in the literature (see p], [8]). These are the ones with the 
property 

d 

2u u >Y,\u ij \, i = l,-.,d. (1.8) 

i=i 

Let D be the set of symmetric matrices satisfying (11.8P and such that 
trace u = 1. The author heard some doubts that, say the results of [3] are 
applicable to equations whose variable coefficients of second order deriva- 
tives form matrices of class D. The point is that the equations in [3] are 
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assumed to have the structure associated with (|1 .3f) with Lipschitz continu- 
ous <Tfc. A naive way fails to take ei,...,ed as the standard basis vectors in 
M. d and write a representation like (jl.4p for a -D-valued function u(x) as (see 

EU) 

u{x) = - \u»{x)\)ei^ + (l/2)(^'(x)) + ( ei + e,)(e t + ej)* 

#3 

+(l/2)(«y(x))-(e i -e i )(e i -e i )*], 

where we used the notation a ± = (l/2)(|a| ± a). The reason for the failure 
is that no smoothness assumptions on u(x) can guarantee that [(u i3 ) + ] 1 ' 2 is 
Lipschitz continuous for i ^ j. One needs a nontrivial structural assumption 
for that. 

Nevertheless, in [7] for d = 2 the author gave explicit formulas for rep- 
resenting twice continuously differentiable D-valued functions in the form 
(jl.7p with Lipschitz continuous p k (u(x)). The result of the present article 
shows that such representation exists for any d. In addition, it turns out 
that the set of diagonally dominant matrices can be replaced with any set 
which is a polyhedron in the set of d x d matrices. By the way, observe that 
(jl.8p can be equivalently written as 

d 

2u" > ^e ij u ij , i = 1, ...,d,e ij = ±1. 

3=1 

In D we also have tracer = 1. Therefore, the bounded set D is described 
by means of finitely many linear equalities and inequalities, and hence D is 
a polyhedron in the space of d x d matrices. Speaking about the case that 
d = 2, it is also worth noting that in [2] an efficient algorithm is introduced 
for approximating arbitrary 2x2 nonnegative matrices with matrices of the 
form ^kPkCkCt, where £l 2 . 

Finally, we reiterate that representation (II. 7h leads to (11.3|) and the latter 
means that we have the following factorization: 

u = vv*, where v tk = akj tk - 

Starting from this point we forget about matrices and work with functions 
having values in a polyhedron. Our main results are presented in Section [21 
Theorem 12. II is proved in Section [2] and Theorem 12.21 is proved in Section HJ 
Section [3] contains an investigation of an auxiliary function some additional 
information on which is provided in Section [5j 

The author discussed the article with Hongjie Dong whose comments are 
greatly appreciated. 

2. Main results 

Let P be a closed bounded convex polyhedron in W d with distinct vertices 
ai, a n , where n > 2. Let d\ > 1 be an integer. 
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Theorem 2.1. Lei 6e a P -valued functions defined on M. dl . Assume 
that the first and second order derivatives of u are bounded and continuous 
on M. dl . Then there exist real valued nonnegative functions ui(y), u n (y) 
such that 

^•Ufc(y) = 1, u(y) = ^2u k {y)a k , (2.1) 

k k 
1/2 /] 

and u k are Lipschitz continuous on R 1 with a constant which depends only 
on P and swp{\u^^(y)\ : \rj\ = l,y G R dl }. 

Clearly the following assumption which we keep throughout the paper 
does not restrict generality: 

a\ = 0, Span (02, a n ) = R d . 

To prove Theorem 12. II we need the following result. For £ G R rf \ {0} and 
x G P denote by ei(a:,£) the distance from x to 5P along the ray x + 
t > 0. Introduce, P° as the interior of P. 

Denote by $ the set of d — 1-dimensional faces of P and for T € $ and 
x G P° introduce dr(x) as the distance from x to T. Also let rip be a unit 
normal vector to T. 

Theorem 2.2. On P there exist Lipschitz continuous nonnegative functions 
Pi(x), ...,p n (x) which are infinitely differentiable in P° and such that 

(i) p k > in P°; 

(ii) in P we have 

y^Pkjx) = 1, x = y^ j p k {x)a k ; 

k k 

(Hi) for any £ G R rf and x G P° we have 

^£^< Nm JS^A, (2 . 2) 

where N is a finite constant depending only on P. 

Proof of Theorem 12.11 Take a point xq G P° and for t G (0, 1) set 

ut{y) = tu(y) + (1 — t)xQ. Then ut takes values in P°. Assume that for 
P°- valued functions the statement of Theorem 12.11 is true. Then, for each 
t G (0, 1) there exist real valued nonnegative functions u t ±(y), ...,u tn (y) such 
that 

5^wtjfe(y) = i, My) = ^2 u tk(y)ak, 

k k 

and u^ 2 are Lipschitz continuous on R dl with a constant independent of t. 
By the Arzela-Ascoli theorem it follows that there exists a sequence t n f 1 

1 /2 

such that ut n k(y) converge to some functions itjt(y) for each y and u k are 
Lipschitz continuous. Obviously, these are the functions which we need. We 
see that without losing generality we may assume that u(y) G P° for all y. 
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We will be using the well-known fact that if we have a nonnegative twice 
continuously differentiable function f(y) given on M. dl and having bounded 
second-order derivatives, then for any y G M. dl 

|V/(y)| 2 < Af(y) S np{\f {v)M (z) : \r,\ = l,z G R dl }- 

Now take pk from Theorem l2.2l and set Uk(y) = pk(u(y)). Then the equa- 
tions (|2,ip obviously hold. Since pk are positive and infinitely differentiable 
in P°, Uk are positive and continuously differentiable in Mr 1 . Therefore, to 

1/2 

estimate the Lipschitz constant of u k it suffices to estimate its first order 
directional derivatives. 

Fix a y, r\ G M dl with \r)\ = 1 and set x = u(y), £ = u^(y). Then by (|2.2p 

o\( l/a ^ f M ^ , r l(™r,OI /r) q\ 

2IK )(„fe)l -j^-fgf ^ • < 2 - 3 ' 

Next, take a face r £ f and let it be given as {x : (rer, x) = b}, where b is 
a constant. By multiplying nr and b by —1 if needed we may assume that 

(nr,w) >b Vu> G P. 

Then f(z) := (nr, u(z)) — b is a nonnegative twice continuously differentiable 
function on M. dl . By the above 

l(nr,0! 2 = \f {v )(y)\ 2 < Nf(y) = N\(n r ,x) - b\ = Nd r (x), 

where 

N = 4 sup \u iv)iri) (z)\. 
This and (|2.3[) bring the proof of the theorem to an end. 

3. AN AUXILIARY FUNCTION 

For x G P° define 

n n n 

U(x) = max{^lnpi : p { > 0, ^Pi = 1, = 

i=l i=l i=l 

Obviously, U < and for each x G P° there exists p\, ...,p n achieving the 
maximum. 

Lemma 3.1. (i) The function U is strictly concave and therefore continuous 
in P°. 

(ii) For each x G P° there exists a unique set pi(x), ...,p n (x) > such 
that 

n n 

),Pi{x) = 1, S ^p i {x)a i = x, U{x) = ylnpijx). 

i=l i=l i 

(Hi) The functions pi (x), ...,p n (x) are continuous in P° . 
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Proof, (i) Take x, y 6 P°, t,s G (0, 1), such that t + s = 1, and let pi, ...,p n 
and gi, q n be some sets achieving the maximums for x and y respectively. 
Then for = ipj + sg^ we have 

i i 

Hence, 

£/" (ix + sy) > m (*Pi + s Qi) > t y~] lnpj + s ^ In g^ 

i i i 

= «7(x) + stffc/), 

where the second inequality is strict if pi ^ qi for at least one i. This is 
certainly the case if x ^ y, which proves (i). Another case would appear if 
x = y and we assumed that there are two different sets pi, ...,p n and q±, q n 
achieving U(x). But then the above computations would lead to a wrong 
conclusion that U(x) > U(x). This proves (ii). 

Finally (hi) follows from the continuity of U(x) and assertion (ii). The 
lemma is proved. 

Lemma 3.2. The function U is continuously differentiable in P° and for 
any x € P° and £ G M. d , which is represented as 

£ = ^2 Qk{x - a k ) (3.1) 
k 

with some numbers q k , we have 



In particular, as ^ = x — a k , 

1 



Pk(x) 



n-U ix _ ak) {x). (3.3) 



Proof. Fix an xo € P° and let A € M d be such that the graph of the 
function (A, x — Xq) + U(xq) is a supporting plane for the graph of U(x) at 
(xo, U (xo)). Set b := J2 k Ik and write 

xo + t£ = ^[(1 + bt)p k (x ) - tq k ]a k . 

k 

For sufficiently small t we have (1 + bt)p k (xo) — tq k > and 

+ bt)p k (x )-tq k ] = 1. 

fc 

It follows that for small i 

i(A, + 17 (x ) > *7(x + t£) > ln [(! + ^)pfc(xo) " tq k ] 
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with equalities instead of the inequalities for t = 0. By differentiating at 
t = the extreme terms we find 

If there is another vector such that the graph of the function (/j, x — 

xq) + U(xq) is a supporting plane for the graph of U(x) at (xq, U(xo)), then 
the above formula implies that A — fi _L ^. This holds for any ^ admitting 
representation (|3.ip with xq in place of x. Since 

Span{(x - a\) - (x - a 2 ), (x - ai) - (x - a n )} = R d , (3.4) 

any £ has the said property, and hence A = /i. 

Thus, for each point xo E P° there is only one supporting plane at 
(xq, U(xq)). This and the concavity of U implies that U is continuously 
differentiable, A = VU(xo), and the lemma is proved. 

Corollary 3.3. Take any representation 

x = '^q k a k with y~]q k = l. 



Then 

Qk 



E 



, Pk{x) 



n. 



Indeed, it suffices to observe that £ = in (|3.2p . 
Lemma 3.4. The functions U , p k are infinitely differentiable in P° . 
Proof. Denote X{x) = VU(x). Then 



Pk(x) 



and X(x) satisfies 



where 



n- (x -a k ,X(x))' 
F(X(x),x)=0, 



F(X,x) = (F i (X,x),i = l,...,d), F i (X,x) = V - 4). 

n- {x - a k ,X) 

We have 

~ X - (.-(,-«.. A))' ^ - - <»• 

k 

By (|3.4p there is no nonzero vectors rj that are orthogonal to all x — a k . It 
follows that the matrix with the entries ^|jF*(A,x) is nondegenerate, X(x) 
is infinitely differentiable by the implicit function theorem and the lemma 
is proved. 
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Lemma 3.5. Let x € P°, £ € R d . Then 

tfftxofr) = - 2^ p2 (a .) ■ ( 3 - 5 ) 

Furthermore, if t; = Y^kQk(ak ~~ x ) f or some numbers q k) then 



W^ir'^^' (3 " 6) 



Finally, \Pk(a k ~x)( x )\ < (n + 1) 1 / 2 /or any fc = l,...,n. 
Proof. By differentiating f)3.3|) we find 

-pTT^y = -*)(€)(*) + U {0 (x). 

By multiplying this equality by q k and summing up with respect to k we get 
(|3.6p provided that £ = J2k ?k( a fc — %)• Differentiating Y2kPk( x ) a k = x and 
T,kPk( x ) = 1 yields 

fc A; fc 

which allows us to use (|3.6p with ^ = an d obtain f)3.5j) . 
Next, the right-hand side of (|3.6p equals 



Its square by Holder's inequality is less than 

e|e(%^- 2 ^)%)+^))- 

We recall (j3.5[) and observe that 

E^) = ' E^ 1 * ^) = (E^) 2 ^ 



Then we find that 



2 



^(k)^ + 1 )l%)(dE§ 

which is equivalent to (|3.7p . 

The last assertion of the lemma is obtained by taking £ = a k — x in (|3.7p . 
The lemma is proved. 
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Theorem 3.6. For £ 7^ in P° we have 

U m )W > -(" + 4w2 ) d 2 (g>0 ^ d 2 (3;> _ - 

In particular, for any x in P° we have 

\PHd x )\ <N _ ISI 



p k {x) d(x, f) A d(x, -£) ' 

w/iere N = (n + An 2 ) 1 / 2 . 

Proof. Without losing generality we assume that |£| = 1, take x € P°, 
and set 

i/ = s + d(a;,Of, I/(t) = (1 - t)x + ty, te {-e,l), 
where e > is to be chosen later. Certainly there is a representation 

Therefore, for sufficiently small e and all i € (— e, 1) we have 

y(t) = ^2 OfcCC 1 ~ ^jPfcfc) + C 1 - OPfcfa) + % > 0, 

k 

J2((l-t)Pk(x)+tq k ) = 1. 

k 

By definition, 

with equality for t = 0. Therefore, the second derivatives in i at t = of 
the extreme terms are linked by a similar inequality, that is 

1 (x, W*) > - £ ^ > -n - £ 
In like manner for z = x — d(x, — £)£ we find 

Z = Y rkCLk > Tk - °> rfc = 



A' 



d 2 (x, -^)U im) (x) >-n~Y^ 



p 2 k (x) 



k 

However, for some a, (3 > such that ay + (5z = x and a + (3 = 1 we have 



^2(aq k + /3r k )a k , ^(ctgfc + /3r fc ) = 1. 



By Corollary E31 

ag fc + /?r fe 



PfcO*0 



n. 
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It follows that 



and hence at least one of 



k > k 

is less than 4n 2 . This yields the result and the theorem is proved. 

Now we are going to get prepared to estimating the Lipschitz constants 
of Pfe's. Recall that a% = and let Pi be the polyhedron with vertices 
02, ■■■,a n , let U n -i(x) be the function U defined relative to Pi, and let P° 
be the relative interior of Pi . 

Lemma 3.7. Let x G P° and let A > 1 be such that Xx G Pf . TTien 

U n -i ( Ax) - n In A + ln( A - 1) < U(x) ; (3.8) 

(ii) we have an equality in (|3.8p instead of the inequality if we take A = 
\{x) := (l-pi(x))- 1 . 

Proof. Let p2, ■■■,p n be the set that achieves U n -i(Xx). Then 

n n 

y^^pi = i) E^ a * = ( 3 - 9 ) 

i=2 i=2 

Therefore, for p 4 := \~ 1 pi, i = 2, n, and p\ := 1 — P2 — ■■■ — p n = 1 — A -1 
we have p\ > 0, since A > 1, and 

n n n 

i=l i=l i=2 

By adding that 

U n -x{\x) = ^\npi = (n- l)lnA + ^lnpi - ln(l - A" 1 ). (3.10) 

i=2 i=l 



we certainly obtain ()3.8j) . 

To prove assertion (ii) observe that for A = A(x), pi := Pi(x), and p~j := 
Apj(x) we have (|3.9p and p~j > 0. It follows that Ax G P° and the first 
equality sign in (|3.1U|) should be replaced with >. By combining this with 
we get what we need. The lemma is proved. 



Corollary 3.8. For x G P° we have A(x)x G P° and the set \{x)p2{x) , . . . , 
X(x)p n (x) achieves U n -i(X(x)x), so that if for y G Pf we denote byp%{y),..., 
Pn(y) the set that achieves U n -\{y), then for x G P° and k > 2 we have 

p k {x) = (1 - pi(x))p k { X ). (3.11) 

1 -Piix) 
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Theorem 3.9. The functions pi(x), ...,p n (x) are Lipschitz continuous in 
P° and, therefore, admit extensions to Lipschitz continuous functions in P. 

Proof. We will be using the induction on n. If n = 2 and say P = [0, 1], 
a\ = 0,a2 = 1, then pi(x) = 1 — x,p2(x) = x, and our assertion is true 
indeed. 

Assume that our assertion is proved for all polyhedra with n — 1 vertices. 
Then the functions pt introduced in Corollary 13.81 are Lipschitz continuous 
in P°. Since |pi(a.)| < (n + l) 1 ^ 2 , from (|3.1ip we have that for k > 2 and 
x G P° 

\Pk(x)i x )\ < N + (! -Pi(x))lim^ 1 |p fc (- — ■ -) -p k { 



ej.o 1— pi(x + ex) l—pi(x)' 



< N + N(l -pi(x))lime 



_ 1 , x + ex 



e[o '1 — pi(x + ex) 1— pi(x)' 

= AT + iVx 1+ ; {) , \ \ <N + N-^—. 

l-pi(x) 1 l-pi(x) 

Thus, for any e > 0, are bounded as long as |x| > e and k > 2. 

Above we also used that P\< x ){x) is bounded. 

Generally, Pk(x-a )( x ) are bounded as long as \x — aj\ > e. In particular, 
Pi(x-a,j){x) are bounded for \x — a,j\ > e. We now claim that there exists 
an e > and No such that, for any unit £ G M. d and x £ P one can find 
numbers 771,...,% such that 

f = n k{x-a k ), \r) k \ < N . 

k:\x- afe|>e 

Indeed, if e is small enough the restriction of summation may exclude only 
one term with k such that \x — a k \ < e. Still the remaining set {x — aj,j ^ k} 
would be close to {aj~ — aj,j 7^ k} a subset of which forms a basis in M. d . On 
the other hand, if there is nothing to exclude, our claim follows from (|3.4[) . 

This proves that |pi(£)(x)| is bounded for x € P°, |£| = 1. Of course, the 
same holds for other |pfc(^)(x)| and the theorem is proved. 

4. Proof of Theorem 12.21 

First we introduce a few new objects. Let an integer r G [2, d] and let 
Ti, ...,T r € $ be such that nr 1 , ■■■,n-p r are linearly independent. Then 

l n r r - nspan^^.^n^jnrvl > 0, 

where Il£ is the orthogonal projection operator on a subspace C G R rf . Since 
there are only finitely many such r and T\, T r G <&, we see that there is a 
constant k > 1 such that we always have 

l"r r -IIspaii(nr 1 ,...,nr r _ 1 ) n rrl > ^ ■ 
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For a Ci > define recursively 

r-l 
i=l 

Obviously, ( r > Cr-i and Cr is a linear function of Ci so that we can choose 
and fix a (j > such that 

7r : = 2 r (r < 1/2, r = l,...,d+l. 

Set 

e r = K _1 7 r (< 7 r ). 
Now fix x and drop it in some notation. Take a V £ $ for which d(x, £) A 
d(x, — £) equals the distance from x to T along the line x + i^lCl" 1 , i £ 
Then 

^,OA^,- f ) = |«|jg|. 

Denote 5 = | Vp^ I — 1 - There are two cases. 
Case 1 

|("r,£)| ^ £ i 5 bfc(e)l- 

Case 2 

|(n r ,OI < £ i 8 \Pk(g)\- 
In the first case additionally assume that Pk{x) < dr(^)- Then by Theo- 
rem [376] 

Kol pI /2 = ri /2 |(nr,g)| |(n r ,Q| 

pj/a " ia d(*,£)Ad(*,-e) d r " 4/2 • 

On the other hand, if Pk(%) > dr^), then 

|Pfe(g)l < (e ^-i l("r,Q| < , ^-i l(n r ,Q| 
~W ^ ( £ld > — 172- ^ ( £ld > ,1/2 • 
ft ft d r 

Thus, in Case 1 we have 

W < WM )-H^|M, 

ft "r 
which proves (|2.2p since 5 _1 is a bounded function. 

In the rest of the proof we concentrate on Case 2. We will be using a 
recursive procedure. Denote £1 = £, Ti = T, ni = nr 1 , and introduce 

£2 = 6- ni(rai>£0- 

Observe that (Case 2) 

_ ft(£i)l = bfc(m)l • < £ _1 |(«l,£l)| < £ l|ft(6)l ^ 7l|jPfc(eol, 

(l-7i)b*(€i)l < |Pfc(&)l < (H-7i)IP*(€i)l- 
In particular, Pfc(£ 2 ) ^ and £2 7^ 0. Also in Case 2 we have £2 -L "-1 and 

16 - £i| = I < £ i^bfc(6)l - 7i^b*(5i)l- 



SMOOTH FUNCTIONS WITH VALUES IN POLYHEDRA 15 

It follows that for an integer r > 2 we have I\ 6 i = 1, ...,r — 1, and 
vectors £j ^ 0, i = l,...,r, such that for nj being the normal vectors to I\ 
we have 

(i) m,...,n r _i are linearly independent; 

(ii) for 1 < j < i < r we have £j _L nj; 

(iii) for 2 < i < r we have 

(l-7i-l)bk(6_i)l < bfc«i)l < (l + 7i-i)b*(€i-i)l; 

(iv) for 2 < i < r we have |& - < 7j_i<5|p fc ( f |; 

(v) for 1 < i < r - 1 we have |(nj,&)| < £i£|p fe (£.)|: 

(vi) for 1 < i < r— 1 the face Tj G is the one for which d(x, £$) Ad(x, — 
equals the distance from x to Tj along the line x + t € K. 

In light of (i) we have r — 1 < d. Also observe that by virtue of (iii) and 
(iv) for 2 < i < r (recall that < 1/2) 

r-l 

\p m _ l} \ < \p k{& )\ II C 1 - ^ ^ 2M+ W)I> 

j=i-l 

r r 

|£r - Cll < E 1^ " ^-i I ^ ^bfcttr)! 5>i-i2 r - i+1 = er5b fc (^)|/2. (4.2) 

i=2 i=2 

Now introduce T r as the face of P for which d(x,£ r ) A <i(x, — £ r ) equals 
the distance from x to T r along the line x + ^r|Cr| _1 ; t £ R. Set n T = nr r 
and first suppose that 

|(nr,Cr)| > £rS\Pk(£r)\- ( 43 ) 

Then as in Case 1 

^<iY M )-lfe|M. (4.4) 

Pk d T r 

Here by (|4.ip the left-hand side dominates 

9 - r bfc(|0l 

Z 1/2 • 

To estimate the right-hand side of (|4.4p use (|4.2p and (|4.3p to get 

|(n r ,Cr)-(n r ,Ci)| < 1^-61 < (1/2)1(^,^)1- 

Hence, 

|("r,6)l <2|(n r ,a)| 
and going back to (|4.4p we obtain 



1/2 - ,1/2 ' 

Pk d r r 

which proves (|2.2p . 

In the situation that (14.31) is violated introduce 



Cr+1 Cr fri^riCr 
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where f r = h r /\h r \ 2 and 

h r = Tl r — ^-Span(ni,...,n r -i) n r- 

Observe that h r ^ 0. Otherwise, n r would lie in Span (m, n r _i), £ r 
would be orthogonal also to n r and the line x + t£r|£r| _1 > t 6 I, would 
have never met r r . In particular, property (i) holds with r + 1 in place 
of r. By the definition of k we have \h r \ > Then |/ r | < « and since 

\(n r ,Cr)\ < e r 5\p kl g r )\, we have 

implying that Pfc(£ r+1 ) 7^ 0, £ r +i 7^ and (iii) holds with r + 1 in place of r. 
Also notice that, for j < r — 1, we have £ r ,h r ,f r _L raj, which implies that 
£ r +i _L n,-. Furthermore, 

(/i r ,n r ) = (h r ,n r - Proj Span(nij . )nr l) n r ) = |/i r | 2 , (/ r ,n r ) = 1, 

(£ r+ l,n r ) = (£ r ,ra r ) - (f r ,n r ){n r ,£ r ) = 0, 

so that f>+i _L rij for all j < r and (ii) holds with r + 1 in place of r. 
Properties (v) and (vi) hold with r + 1 in place of r by the assumption and 
construction. 
Finally, 

|fr+l -£r| = |/r| " |(»r,£r)| < K£ r 5|p fc ( 5r ) | = 7r£bfc(f r ) |, 

so that (iv) holds with r + 1 in place of r. 

Thus, if (|4.3p is violated, we can find objects 1^, £j 7^ having the proper- 
ties (i)-(vi) with r + 1 in place of r. This recursive process will stop at least 
when r reaches d + 1, just because property (i) will prevent us from finding 
n r+ i, which implies that at least at this moment (j4.3l) should be satisfied. 
This proves the theorem. 

5. Additional information 

Remark 5.1. One can estimate p k {x) from below for x 6 P°. It turns out 
that 

, \ ^ d(x,x- a k ) 

Pk{x) > —77 r- — ; r, (5.1 

nd(x, x — ak) + n\x — ak\ 

where nd(x, x — a k ) + n\x — a k \ is obviously bounded away from zero. This 
and the fact that |V»&(:c)| is bounded, actually, show that Pk(x) behaves 
like d(x, x — a k ). 

Indeed, take any £ G M. d \ {0} and observe that y := x + £d(x, ^ P 

can be written as 

fc A; 

Then _ 

£d(x, = y-x = ^2(pk(x) ~ Qk)(x ~ a k ) 

k 
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and by fl32J 

u (Q {xw Xi o/iei = E -^p = E ^ - -n, 

so that, for any £ / 

If £ = x — a/;, this and (|3.3p imply that 

1 nix — aJ 

_ _ = n - f7 (a ._ a )(x) < n + — -, 

Pk{x) d(x,x-a k ) 



which is equivalent to (|5.ip . 

Remark 5.2. If £ = (^(a^ — x) and q% = then 

V- Jl_ = V- (Pfc(g)(^)) 2 + (Pfc(fl(aQ -9fc) 2 (5 2) 



Indeed, write 



Zv„2^ ^2^ ~^Z^ „2/^N „,2, 



-^p^) k pi» pLu^) v 

and observe that the middle term on the right is zero due to (|3.5p and (|3.6 



Remark 5.3. One can improve the estimate of \p k { x -a k )\ from Lemma 1331 
It turns out that 

1 > Pk(x-a k ) + 1 - Pk = OikPh, ( 5 - 3 ) 

where 

(Pk(x-a k ) + 1 - PA:) 2 V- (Pi(x-a k ) ~ Pif 

Ct k = 2 1" 2 • 

Pi 

In particular, p fc > p k ( x -a k ) > Pk ~ 1 - 

Indeed, we may concentrate on proving (|5.3I) only for k = 1 in which case 
we apply (|5.2p with £ = x — ai . Since 

x-a\ = (pi- l)(ai - x) +^2(02 - x) + ... + p n (a n - x) 

one can take q\ = p\ — 1, q2 = P2, Qn = Pn- Then by (|3.5p and (|5.2p (recall 
that a\ = 0) 

, , (Pi-l) 2 r , 
Pi 

On the other hand, by differentiating (|3.3p we find 

tt , tt tt , 1 
—p- ~ u (x)(x) + t>(x) - f(x)(x) + n - — • 
Pi Pi 

Hence, 

1 , (Pi ~ I) 2 Pi(x) , 1 , 
n — 1 H ^ = ^ — h n h ai, 

Pi Pi Pi 
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-2pi + 1 = -pi( x ) -pi + oiip\, 



and the equality in (|5.3[) follows. The first inequality follows from the fact 
that a x p\ > (pi (a .) + 1 - pi) 2 . 

Remark 5.4. On can combine Remark 15. II and the fact that for any j, k 

\Pk(x-a/)(x)\ <(n+ l) 1/2 p fc (x)/p i (x). 

in order to investigate the behavior of Vpi(x) as x approaches dP. The 
above mentioned fact follows from (|3.7p when £ = x — a,j. 

Remark 5.5. The functions p^ have a peculiar symmetry. It turns out that 
for all j, k 

Pk{x-a 3 ) 1_ = Pj(x-a k ) J_ 
Pfc Pi ^ Pfe' 



Indeed, differentiating (|3.3p easily yields 

P * ( T 3 ° + ^ " " = U (*-**)(*-°s)> ( 5 - 4 ) 
Pfc P? 

where the right-hand side is symmetric with respect to j, k. One can com- 
bine (|5.4[) with Remarks 15.11 and 15.41 to obtain some information about the 
behavior of the second-order derivatives of U near the boundary of P. For 
instance, \U^ ak ^ x _ a .)\ < n + 2(n + l) 1/2 (pfcPj) _1 . 

Remark 5.6. Lemma 13.21 allows one to obtain a precise information about 
the behavior of the first-order derivatives of U near the vertices of P. Indeed, 
(|3.3p shows that 

lim Ui x _ a )(x) = n - 1. 

Remark 5.7. If in the situation of Lemma 13.71 we have X(x)x G P° then the 
derivative with respect to A of the left-hand side of (13. 8h is zero at A = A(x). 
In this case, by substituting A = A(x) into (|3.8p and differentiating with 
respect to x we find 

\(x)U n _ 1>x i(\(x)x) = U x i(x), i = l,...,d. 

In particular, the gradients of U n -\ and U are proportional at corresponding 
points. 
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